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We show that the linearized phase space flow around a discrete breather solution is not capable of 
generating persistent energy flow away from the breather even in the case of instabilities of extended 
states. This holds both for the classical and quantized description of the flow. The main reason for 
that is the parametric driving the breather provides to the flow. Corresponding scaling arguments 
are derived for both classical and quantum cases. Numerical simulations of the classical flow support 
our findings. 
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I. INTRODUCTION 



Discrete breathers (DB) have been intensively studied in the past decade. They are known to be generic solutions 
of the dynamics of nonlinear spatially discrete translationally invariant Hamiltonian systems. DBs are time-periodic 
and spatially localized excitations and belong to one-parameter families of solutions of the underlying equations of 



| DBs have been observed in various experimental situations ranging from Joscphson junction ladders^, coupled 
nonlinear optical waveguides^ and driven micromechanical cantilever arrays^ to layered antiferromagnets^ and high- 
T c superconductors^, surface and bulk lattice vibrations of solids 6 and Bose-Einstein condensates loaded on optical 
lattices^. DBs are also predicted to exist in the dynamics of dusty plasma crystals^. The characteristic spatial scales 
range from micrometers down to Angstroems. Especially in Josephson junction ladders DB excitations have been 
t-H , studied very extensively, including their interaction with the modes of the lattice part which is not excited by the 
£> ■ DBiS. Thus the issue of stability of DB states, of their interaction with lattice modes, and possible mechanisms of 
| radiation of energy by DBs due to this interaction becomes an important and timely issue. 

The fact that DB solutions are generic for nonlinear Hamiltonian lattices implies that in the majority of cases the 
underlying Hamiltonian equations of motion are not intcgrablc. This puts limitations and complications on the study 
of perturbed DB states. If the perturbation is considered to be of small amplitude, its evolution can be described using 
' a linearization of the phase space flow around the DB. The resulting linear coupled ordinary differential equations 
, (ODE) with time-periodic coefficients can be studied within the framework of Floquet theorjiii. 

■ Taking into account higher order terms in the phase space flow will ultimately lead to more complicated nonin- 
tegrable equations, which can be studied only approximativelj*i£. Several issues are at stake when discussing the 

H evolution of perturbed breathers. One can consider localized or extended perturbations on one hand. On the other 
hand there are differences in the way breathers react to such perturbations depending on the amplitude of the latter. 

Let us first discuss the case of a linearized phase space flow around a DB. Formally the obtained Floquet equations 
decouple the dynamics of the DB (which is assumed to be given) from the evolution of the perturbation. Perturbations 
which grow in time will then invalidate the abovementioned linearization. Perturbations which decay in time do not 

■ contradict the linearization, but in fact the energy stored in the initial perturbation cannot simply disappear if we 
consider Hamiltonian dynamics. Consequently there is a subtle way this energy will have to be transferred to the 
DB, which again is beyond the linearization frame. Marginally stable perturbations, which neither grow nor decay 
seemingly do not violate the assumed linearization. Nevertheless it has been shown that for extended perturbations 
such a case may be accompanied by a nonzero energy flux emitted out of or into the breather core—. Again it 
would violate the linearization frame. Even though it may do so, such predicted radiation scenaria are confirmed in 
numerical simulations, underpinning the use of the linearization picture. 

Additional sources of radiation can appear when taking into account nonlinear corrections to the phase space flow 
of a DB perturbation, even if the linearized case did not provide with such sources. For instance a marginally stable 
localized linearized perturbation will yield an energy radiation due to the appearance of new frequency combinations 
and resonances with the spectrum of small amplitude plane wavesii 

Here we will be concerned with a particular case within the linearized phase space flow frame, which corresponds 
to the abovementioned extended perturbations yielding a nonzero energy flux out of the breather core. The question 
we want to pose is whether this energy flux can be sustained if the perturbation we choose is local in space. The 
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perturbation will have some overlap with the extended ones. So there will be some radiation, but at the same time 
the initial localized perturbation will simply disperse away from the breather lowering its amplitude. The question 
then is whether these two counteracting processes balance each other or not. The question is of relevance also in 
connection with recently discussed radiation mechanisms of strongly excited quantum breathers^. We will provide 
with answers for both cases. 



II. THE CASE OF CLASSICAL BREATHERS 



For the sake of simplicity we consider first a one-dimensional lattice with the equations of motion 

xi + V'{xi) + W'(xi - xi-!) - W'(x l+ i -xi) = 0, (1) 
which corresponds to the Hamiltonian 

'I 



2 ±f + V{X l ) + W(x t - xi-!) 



(2) 



The index I denotes the lattice site, and can run over a finite or infinite lattice. Extensions to higher lattice dimensions 
are straightforward. A local minimum energy state xi — X\ — is provided by V(0) = W(0) = V'(0) = W'(0) = 
and V"(0),W"(0) > 0. Small amplitude excitations can be obtained by linearizing Q with the ansatz 

a>j(t)~e i < &, «*-« , > (3) 
which results in the dispersion relation for plane waves 

^ = y"(0)+4W r "(0)sin 2 | . (4) 

For varieties of anharmonic potentials V, W it is well knowni that the equations QJ allow for families of discrete 
breather solutions of the type 

x t (t) = xi(t + T b ) , xi^oc -> . (5) 

Here the breather frequency fl b = 2ix jT b is a tunable parameter which satisfies the nonresonance condition 

mfl b 7^ oj q . (6) 

In addition given a breather family we can generate new families by discrete translations I — > I + Zo- 

In the next step we consider small perturbations ti around a given breather solution for a given value of its frequency 
fl b . We insert the ansatz xi(t) — xi(t) + ej(t) into the equations of motion (|TJ and linearize it with respect to the 
perturbation: 

h + V"{ Xl (t))ei + W"(xi(t) - xi-x(t))(ei - e,_i) - W"{x l+1 {t) - Xi(t))(e l+1 -ei) = 0. (7) 

Integration of these Floquet equations over one breather period T b maps the phase space (ez,ez) onto itself and is 
equivalently described by a Floquet matrix F. This matrix is symplectic and can be obtained e.g. numerically. Its 
eigenvalues A and eigenvectors y\ describe the linear stability properties of the DB and the scattering of plane waves 
by the DB as well. If |A| = 1 the corresponding eigenstate y\ is marginally stable - it neither grows nor decays in time. 
Eigenvectors may be localized or extended. Since the breather is localized, the eigenvalue spectrum of the localized 
eigenstates is discrete, while the eigenvalue spectrum of delocalized eigenstates is continuous for an infinite lattice. 
Indeed, in that case delocalized eigenvectors far from the breather asymptotically take the form @. Consequently 
their eigenvalues are given by 

A, = e iUJ " T " . (8) 

Note that this is true only for an infinite lattice. In the following we will remind the reader about finite size corrections 
to this picture. 

Let us now consider a lattice with a particular breather solution such that 

L0 q ± uo q i = mVt b , q ^ q' , (9) 
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so that X q — \ q i. This twofold degeneracy is the origin of inelastic multichannel scattering 13 performed at these 
wavenumbers q and q' . Note, that such a situation corresponds to the parametric resonance in the system (JJJ, where 
the breather acts as a parametric driver (exponentially localized in space). Thus an important question arises: can 
one pump energy into the system in the regime of a parametric resonance, provided that parametric driving is local 
in space? The answer is no, as long as we deal with an infinite system size. Here however we note, that for a finite 
lattice the outcome will be close to the above statements, but not identical. Indeed, as shown by Marin and Aubrj*i&, 
the degeneracy of the corresponding Floquet eigenvalues is lifted in a finite lattice, leading to their departure from the 
unit circle, so that |A| ^ 1. However these departures are the smaller the larger the lattice size is. Marin and Aubry 
have shown using a band analysis that ||A| — 1| ~ 1/N, where N is the total number of lattices sites. This result is 
obtained first by noting that the twofold degeneracy is lifted due to the time-periodic localized DB perturbation and 
can be accounted for by estimating the corresponding matrix element. Because both initially degenerate eigenstates 
are delocalized, whereas the DB is localized, the matrix element will be of the order of 1/N. Perturbation theory of 
two degenerate states tells immediately that the degeneracy will be lifted to an amount of 1/N as well. 

It is instructive to revisit here the consequencies. Assume we have a large but finite lattice and such an extended 
slightly unstable eigenstate. Taking the perturbation along this eigenstate, and integrating over one period of the 
breather Tb, this perturbation will grow everywhere in the system, though not very strongly. How can that happen, 
if the breather itself is localized, say exponentially? The influx of energy is provided by the breather, and is confined 
to a finite part of the lattice. The only possibility is that there is a nonzero energy flux in the outer regions of the 
lattice due to some slightly inhomogeneous profile of the Floquet eigenstate. Then the breather simply feeds energy 
in (or out) of a confined small part of the lattice, and this energy then travels along the rest of the lattice. If so, 
the legitimate question arises whether this can be used as a possible source of radiation of waves by a breather, if a 
corresponding generic localized perturbation of the breather is excited. 

To answer this question let us assume that all the elements of a given (extended) Floquet eigenvector are of the 
same order of magnitude y\ ~ A. Then taking an initial perturbation being equal to such an eigenvector the stored 
energy is given by E\ ~ NA 2 . If we choose a growing unstable eigenvector, after one period Tf, the amplitude of 
the perturbation will become (1 + 1/N) A. Consequently the energy after that time increases to N(l + 1/N) 2 A 2 . 
Thus the energy grows during that time by an amount of AE ~ 2A 2 for large N. Assume now that we make a local 
perturbation of the breather given by some vector B with its elements again being at most of order A. Then it can 
be represented as a superposition of the Floquet eigenvectors: 

£~4$>. do) 

A 

Consequently the growth of energy in each unstable eigenstate during one period Tb of the breather amounts to 
AE\ ~ A 2 /N 2 . The total number of unstable eigenstates will be always less than N, so that the total growth of 
energy in the direction of all unstable eigenstates is limited by A 2 /N. Recalling that A is fixed, this growth rate will 
scale to zero for an infinite lattice. Consequently we predict that the above discussed finite lattice instabilities do not 
result in an energy radiation of the breather when perturbed locally. In other words, a local perturbation will start to 
disperse away from the breather, lowering its amplitude in the breather core. At the same time it has some overlap 
with unstable eigenstates which would result in a growth of the same amplitude. However the dispersion acts more 
efficiently, and for large times the amplitude of the perturbation in the breather core will ultimately decay down to 
zero. 

The above argumentation was provided for a one-dimensional system. It is straightforward to generalize it to higher 
lattice dimensions d as well. Assuming now that N represents the linear size of the system, the stored energy in an 
unstable eigenvector becomes E\ ~ N d A 2 . After one period Tf, the energy growth will be given by AE ~ 2N d ~ 1 A 2 . 
For a local perturbation the growth of energy in each unstable eigenstate per Tb is then given by AE\ ~ 2N~ d ~ 1 A 2 . 
Since there are at most N d unstable eigenvectors, we obtain again an upper limit for the total energy growth A 2 /N, 
independent of the dimension of the lattice. 

In order to verify this prediction, we performed high precision numerical simulations. We consider a system with 
V(x) = x 2 /2 — x 3 /3 and W(x) = 0.5Ca; 2 . First we consider a lattice of size N = 40. We identify a breather with 
frequency fib = 0.75 for which the abovementioned finite lattice instabilities take place in the coupling constant region 
C <J 0.093. We compute the eigenvalues and eigenvectors of the corresponding Floquet matrix for two different values 
of the coupling constant: C = 0.092 (breather is linearly stable) and C = 0.093 (breather has finite size instabilities). 
We then perturb the breather in the direction of a stable (at C = 0.092) and a slightly unstable (at C = 0.093) 
extended eigenstates and monitor the energy flux at some distance from the breather 



jl(t) = e z (e;_i - ei +1 ) 



(11) 
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FIG. 1: Energy flux ji(t) (solid gray lines) and integral energy flux Ji(t) (solid black lines, see also insets) in breather dynamics 
with small perturbation for the cases: (a) C = 0.092, perturbation along stable eigenstate; (b) C = 0.093, perturbation along 
unstable eigenstate. The system size is N = 40 (breather is centered at 20th site), energy flux is monitored at the site I = 25. 
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FIG. 2: Eigenvalues \ q of the Floquet matrix F for the cases: (a) C — 0.02, Qt, — 0.75, inset shows the breather profile (central 
part); (b) C = 0.1, Qb = 0.75, inset zooms the region of instabilities. The system size in both cases is 80 sites. 



as well as its integral over the observed time of simulation r 

Ji(t) = - I j t (jt)dt. (12) 
T Jo 

We nicely reproduce the expected two different scenaria: oscillations of energy flux with zero average of its integral 
value in the case of a stable breather; and slow but still exponential growth of the energy flux out of the breather 
region in the case of an unstable breather, see Fig. 

Next we take a large lattice with N — 4000 and identify breather states with f2f, — 0.75 for two different values of 
the coupling constant: C — 0.02 (the breather is linearly stable, see Fig.|2[a)) and C — 0.1 (there exists a considerable 
number of unstable extended eigenstates, see Fig. Efb)). We excite a local perturbation on top of the breather core 
in a random way on 10 adjacent sites with amplitude being of order one at each site. We then monitor the time 
dependence of the energy flux averaged over 20 lattice sites j sum (t) = J^n + •?"'(*) outside the breather core. In both 
cases we observe a decay of the measured radiation with growing time, see Fig- EI F° r large times it follows the power 
law j ~ i -4 (see insets in Fig. to smoothen the oscillations and make the picture more clear in insets we plot the 
flux value averaged over 4 breather periods: j SU rn(t) = 1/(4T},) J t _ 4T jsum{T)d,T), which is explained in the Appendix. 
This clearly demonstrates that the finite lattice instabilities are not capable of sustaining a nonzero energy radiation 
of the breather induced by a local perturbation. 
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FIG. 3: Energy flux jsum(t) for a local initial perturbation for the cases: (a) C = 0.02; (b) C = 0.1. The system size is 
N — 4000 (breather is centered at 2000th site), energy flux is monitored at sites I £ [2020, 2040], perturbation is made on top 
of the 10 adjacent sites in the breather core in a random way. Insets show the averaged over 4 breather periods integral energy 
flux jsum{t) (see the main body text for details). Dashed line in inset in (a) shows the asymptote j a s(t) = le+8/(t/Tj,) 4 . 



III. THE QUANTUM CASE 

The general problem of quantizing J3J) and understanding the correspondence between quantum eigenstates and 
classical discrete breathers is an issue of current research. Since the underlying models are nonintegrable in general, 
no analytic solutions can be obtained. Approximations have to be complemented by numerical studies. However this 
is in general very hard to achieve, since one has to deal with the diagonalization of huge matrices. Indeed, even when 
studying one single oscillator, strictly speaking one has to consider an infinite dimensional Hilbert space. So already 
at the level of a single oscillator justified cutoffs in its Hilbert space have to be introduced. Taking into account many 
coupled oscillators leads to the necessity to reduce the number of states per oscillator. Reliable results in the high 
energy domain of the quantum problem have been obtained so far only for small systems of two or three coupled 
oscillatorsiLiSii2i2£ . These studies together with computations of larger systems at lower energies have confirmed 
that quantum breather states are many phonon bound states^. They correspond to some extent to classical breather 
excitations being able to tunnel along the lattice^. The tunneling rate is expected (and confirmed for small systems) 
to be exponentially small for large numbers of participating phonons. Then it is legitimate to assume that a localized 
excitation of the lattice with an energy corresponding to a large number of phonons will evolve for exponentially long 
times according to its classical DB analogue. During that time however the evolution of small amplitude perturbations 
around the breather can be considered in its full quantum version. In order to proceed we observe that Q) corresponds 
to a time-dependent Hamiltonian 



(13) 



Here 717 is the canonically conjugated momentum to e;. Now we may consider the corresponding quantum Hamiltonian 
operator 



(14) 



The operators satisfy the standard commutation relation [ej, 7r m ] = idi jTn . 

The fact that the Hamiltonian (|14fl is a quadratic form of operators e/ and 7T; is crucial: we may appeal to the 
Ehrenfest theorem^ and conclude that the dynamics of the quantum system l|14|l is in full agreement with the 
dynamics of the corresponding classical system ff 1 31) and therefore no additional quantum effects which may lead to 
breather radiation should appear. To show this we switch to time-dependent Heisenberg operators ef (t) , nf 1 (t) using 
standard relations: 

A H {t) = (fe i ^ dT ^)A{Te- i ^ dT ^), (15) 

where A H (t) is the time-dependent Heisenberg operator corresponding to a time-independent operator A, T and T 
are time ordering and anti-time ordering operators, respectively. The equation of motion for a Heisenberg operator 
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A H (t) reads: 



f)A H (t) 

-i^-gf 1 = [H e (n,i,t),A H (t)}. (16) 



After substitution of the Hamiltonian (fH|) into (jTT)|) it follows that operators ef (t) satisfy the discussed above 
classical equations Q, in which coordinates ei are substituted by the corresponding operators ef 1 . Since these 
equations are linear we may average them with the time-independent wave function iJjq , corresponding to the initial 
state of the system, and get an equation for the expectation value of the coordinate operator 

e,(t) = (^olef(t)IVo) (17) 

which is identical to with expectation value e standing instead of the classical variable e. Therefore all the 
conclusions made in the previous chapter as for time evolution of the classical coordinate ei hold for the expectation 
value ii of the quantum operator ii as well. 

Alternatively, the question of possible radiation of phonons in the quantum case may be studied by considering a 
quantum Floquet problem. Let us split the Hamiltonian (|14|) into a time-averaged and an ac part: 

H(t)=H dc + H ac {t) . (18) 

The ac part H ac is local in space because of the locality of the discrete breather which is the origin of the ac drive. 
We can consider the full orthonormal basis of eigenfunctions <j> v of the dc part 

H d c4>v = tv<Pv (19) 

and expand the full wave function \P which satisfies the time-dependent Schrodinger equation 

ii> = H(t)W (20) 

in that basis: 

¥(*) = ^C^t)^. (21) 

This will lead to a set of coupled first order differential equations for the coefficients C v : 

iC v = e v C v + J2 h ^) C M ( 22 ) 

where the matrix elements 

M*K<^Ac(i)<M (23) 

have been introduced and (...) denotes the scalar product in the space of <f>. We also note that because H is Hermitian 
(in fact real symmetric) 

^l*(*)| 2 = |£|O,(t)| a =0. (24) 

In order to answer the question of radiation, let us start with noticing that (|22|l constitute a Floquet problem 
similar to the classical case. The difference is that the rank of the Floquet matrix is formally speaking infinite even 
for a finite lattice, since the Hilbert space dimension of a single site oscillator is infinite. The extended states in 
(1 1911 can be characterized by the amount of excited one-phonon energies and classified accordingly as many-phonon 
excitations with a given number of phonons participating. This constitutes the main difference to the classical Floquet 
problem J7J), where only pairs of one-phonon excitations appear (because in the classical case we compute frequencies 
instead of energies, and time reversal symmetry provides with two possible signs of the phonon frequency). For the 
quantum case we have an infinite number of one-, two-, three-phonon excitations etc. To estimate the magnitude of 
the departure of quantum Floquet eigenvalues from the unit circle, we need again, as in the classical case, to first 
estimate the matrix elements. Since the DB solution has a main frequency contribution and higher harmonics with 
amplitudes exponentially decaying with increasing order, the main contribution will originate from the main frequency 
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component of the DB, which couples the space of many-phonon states locally, e.g. the ground state with the two- 
phonon states, the two-phonon states with the four-phonon states, the four-phonon states with the six-phonon states 
etc. In other words, we have to consider an infinite matrix with degenerate diagonal elements and nearest neighbour 
interaction elements of the order of 1/N in analogy to the classical case. The eigenvalue spectrum of such a matrix 
will spread around the value of the diagonal elements to the same order 1/N. Thus we conclude that the quantum 
Floquet eigenvalues will depart from the unit circle not farther than 1/N exactly as in the classical case. 

What remains then is to repeat the final argument applied in the classical case. This can be done in full analogy to 
the previous chapter, noting that the observable of the flux ji(t) will be defined through the corresponding operator j 
and the product j(t) = ( \J/ 1 _7 1 \& } . This product will be a quadratic form of the time-dependent coefficients C v (t), which 
completes the above analogy. The conclusion is thus that despite the fact that the quantum Floquet problem involves 
an infinite number of bands, and despite the fact that the norm is conserved in the quantum case, a local perturbation 
around the breather will not lead to a persistent radiation of phonons. The argument that radiation must take place 
because the ground state of the unperturbed system (without DB) is not anymore the ground state of the system 
with a DB must then be misleading. In fact the computation of the exact quantum Floquet eigenvalue problem will 
show that there is always a locally deformed ground state (as well as the excited states) . That deformation is clearly 
not the cause of radiation. The only possible cause - an instability of extended Floquet states - has been excluded by 
the above reasoning. 



IV. CONCLUSIONS 



In this work we excluded a particular mechanism of radiation of perturbed breathers driven via weak finite size 
instabilities of extended states. While the statement is rigorous when treating the whole system classically, we arrive 
at a similar result also when treating the fluctuations quantum mechanically, leaving the breather solution to be a 
given classical one. One way to obtain nonzero radiation is to include higher order terms of the perturbation e which 
together with possible localized Floquet eigenvectors of the linearized phase space flow will provide with a constant 
radiation rate of the breather into the plane wave continuum, both for a classical as well as a quantum treatment of 
the fluctuations. Another path is to quantize the breather itself in the quantum case. Then we can expect breather 
tunneling along the lattice, which will provide with some diffusion of the full breather energy out of the originally 
excited lattice part. Concluding we may say that discrete breathers are surprisingly robust objects. They can radiate 
energy into the continuum of a large lattice only via higher orders of perturbations around them. 



APPENDIX 

An estimation of a wave packet dynamics, resulted from a local perturbation on a lattice, can be made essentially 
in a similar way as it was done for continuous systems^. Let us start with an instructive case of a lattice without 
any breather. We excite a local initial perturbation, say e;(0) = <5;.o- Its representation in the reciprocal lattice space 
with wave number q is given by a ^-independent constant e q = const. Consequently the evolution of the perturbation 
after some time t will be given by 

e,(t)~ Ce^-^dq (25) 
Jo 

with Uq given by the plane waves dispersion relation Rewriting it as 

e <(9i-w,t) dgr= f e iP(i)tdq (26) 
o Jo 

with F(q) = qj — uj q we can estimate the integral by noting that for large values of t only g-values contribute for 
which dF/dq = k — v q is small. Here v q is the group velocity at wave number q. For large t only waves with small 
group velocities contribute, i.e. with wave numbers q close to the edges of the first Brilluen zone q = q = 0, n. It is 
straightforward to show, that all extremum points q — q* of the function F(q) give contribution to e/(i) of the same 
order in small parameter 1/t and do not cancel each other, hence it is enough to consider only a single extremum 
point. Expanding F around its extremum q = q* , the integral can be estimated to be 

e^dq ~ e tF ^ J ^ F "(i')(<i-i') 2 t dq „ J_ e i(f4+«<-«w*) } (27) 
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where £ = (dv q /dq) 1 \ q= q- Using the definition (111(1 we then obtain 

h{t) ~ t" - 5 (28) 
ei+i(i)-ei_i(i) - r°- 5 sin(g + 2e//t)~t- 1 - 5 (29) 
~ i" 2 . (30) 

We tested this prediction numerically and found complete agreement. 

The observed 1 jt 4 dependence for a perturbation on top of the breather can be now explained by noting that the 
breather represents a local violation of the translational invariance. In such a case the abovementioned reciprocal 
lattice representation e q of a local perturbation becomes g-dependent because plane waves are not the true eigenstates 
of the system anymore. It is easy to show that already for a single site time-averaged breather contribution the 
amplitude of the extended eigenstates at the breather site is proportional to sing. This implies that e q ~ sing. 
Following then again the above reasoning, we obtain 

e,(t) ~ sin(g*)e' iF («*>' / e i0 5F " ^)(i-i') 2 t dq (31) 



m ~ 1- 4 (32) 



confirming the numerical results in FigEI 
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